Low-temperature specific heat of real crystals: Possibility of 
leading contribution of optical and short- wavelength acoustical vibrations 
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We point out that the repeatedly reported glass-like properties of crystalline materials are not 
necessarily associated with localized (or quasilocalized) excitations. In real crystals, optical and 
short-wavelength acoustical vibrations remain damped due to defects down to zero temperature. If 
such a damping is frequency-independent, e.g., due to planar defects or charged defects, these optical 
and short-wavelength acoustical vibrations yield a linear-in-T contribution to the low-temperature 
specific heat of the crystal lattices. At low enough temperatures such a contribution will prevail over 
that of the long- wavelength acoustical vibrations (Debye contribution). The crossover between the 
linear and the Debye regime takes place at T* oc y/N, where N is the concentration of the defects 
responsible for the damping. Estimates show that this crossover could be observable. 



It is generally accepted that, at low enough tempera- 
tures, the specific heat of nonmetalic solids follows the 
Debye law: C oc T 3 Loosely speaking, it can be 
said that the Debye law arises because of the thermal 
activation of the long-wavelength acoustical vibrations 
according to the Bose-Einstein distribution. Optical and 
short-wavelength acoustical vibrations do not contribute 
significantly to the low-temperature specific heat because 
their thermal activation is exponentially suppressed. Ex- 
ception are glasses, for which a linear-in-T dependence is 
known since long ago. Such a linear dependence is usually 
ascribed to the existence of the tunneling two-level sys- 
tems postulated by Anderson et al. and, independently, 
Phillips in the early seventies |2|. Up to date, however, 
the true microscopic nature of these two-level systems is 
unclear. A glass-like specific heat is observed in strongly 
disordered crystals, what upholds the common belief that 
in these crystals the above mentioned systems there also 
exist. "Anomalous" low-temperature specific heats were 
also observed in a number of crystals with relatively small 
concentrations of defects, but not in more perfect crystals 
of the same composition — not in the same temperature 
range at least (see e.g. Refs. 0,0, and the references 
therein). Appart from the evident conclusion that these 
anomalies are related to defects, no reliable explanation 
was found. 

In this Letter, we argue that a linear-in-T specific heat 
is a natural low-temperature property of any crystal which 
contains a small concentration of defects. The key points 
are i) that the optical and short-wavelength acoustical 
vibrations of real crystals remain damped down to zero 
temperature due to presence of defects and ii) that the 
contribution to the specific heat of these damped vibra- 
tions may be linear-in-T at low enough temperatures. In 
the view of this, it is evident that the contribution to 
the specific heat associated with these optical and short- 
wavelength acoustical vibrations may prevail over the 
(acoustic) Debye one at the lowest temperatures. The 
rest of the Letter is devoted to a more detailed argumen- 
tation and to estimations. 

Nowadays, the study of quantum effects in dissipative 



systems is a topic of vivid interest. As a result of this ac- 
tivity, the specific thermodynamic features of the damped 
harmonic oscillator arc now well understood (sec, e.g., 
Ref. [U). It is worth mentioning that these features have 
already been proved relevant to explain some phenomena 
peculiar to vortex lattices in superconductors. We refer, 
in particular, to the significant contribution that vortons 
ield in the corresponding low-temperature specific heat 
: vortons are acoustic phonons in what concerns to the 
dispersion law (they are connected to acoustic vibrations 
of the vortex lattice), but optic ones in what concerns to 
their finite damping at k = (resulting from the vortex 
viscosity). It is worth mentioning that in incommensu- 
rate phases a similar situation takes place: phasons in in- 
commensurate phases are analogous to vortons in super- 
conductors. The low-temperature specific heat observed 
in some incommensurate phases has also a peculiar de- 
pendence on temperature 7] . Discussion on the origin of 
such a dependence will be published elsewhere Q- Here 
we shall concentrate on the low-temperature specific heat 
of "ordinary" crystal lattices. 

A detailed discussion on the free energy of a damped 
oscillator can be found e.g. in Ref. 5]. To our pur- 
poses, it is convenient to start by considering the case of 
a frequency-independent damping. The asymptotic low- 
temperature expansion of the free energy then reads 
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where To is the ground state energy, 7 is the damping 
coefficient, and luq is the natural frequency of the oscilla- 
tor 9] . It is worth noticing the strong difference between 
this power-law asymptotic expansion and the exponential 
one that, in accordance with the Bose-Einstein distribu- 
tion, is obtained for an undamped oscillator This is 
connected to the broadening of the density of states that 
the damping provokes (the broadening is such that there 
is no gap above the ground state, see Ref. (lQ|). 

As we have mentioned, all vibrations of a real crystal 
(acoustical and optical ones) are in principle damped due 
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to defects down to zero temperature. Let us assume for 
a while that this damping does not depend on the fre- 
quency as it is considered in Eq. QJ. At first glance, in 
accordance with this expression, it will be necessary to 
take into account absolutely all vibrational modes of the 
system when further computing, for instance, the cor- 
responding low-temperature specific heat. But a basic 
assumption when obtaining Eq. itself is that a cer- 
tain separation between "relevant" and "irrelevant" de- 
grees of freedom was possible: it is necessary that a part 
of the system acts as a reservoir to produce dissipation 
(damping) [|. If Eq. © is used, some sort of "double 
counting" then seems unavoidable. However this is not 
completely true. Notice that, for a given finite tempera- 
ture, there always exist acoustical vibrations with small 
enough frequencies for which the power-law expansion 
makes no sense (TjJ. As a result, these acoustical vibra- 
tions with small wavevectors give rise, in particular, to 
the well known Debye contribution to the specific heat. 
They are just these acoustical vibrations which act as 
the above mentioned reservoir necessary to provide the 
damping (see below). 

The asymptotic expansion Eq. can then be used 
to further estimate the contribution to the specific heat 
due to "optical" vibrations (optical and short-wavelength 
acoustical ones) in the low-temperature limit. Indeed at 
temperatures T <C hA/ks, where A is the characteris- 
tic frequency of the corresponding "optical" branch, it 
is sufficient to retain the lowest-order term in T. Thus, 
the contribution to the low-temperature specific heat of a 
damped '"optical" branch can be roughly estimated as 
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where d is the characteristic interatomic distance. Here 
the damping constant 7 has been assumed to be the 
same for the whole "optical" branch. This "optical" 
contribution to the specific heat will prevail over the 
acoustic one Q at T < T* ~ ( 7 /A) 1 / 2 6, where 6 = 
he 3 / 2 /(fcsA 1 / 2 ^ 2 ), with c being the velocity of sound, 
is of the order of magnitude of the Debye temperature. 

Let us mention that it is not completely consistent to 
single out one "optical" branch and the acoustic reser- 
voir given that there is a "viscous" coupling between 
these "optical" vibrations due to the damping (described 
by non-diagonal terms in the corresponding dissipative 
function). However, in order to reveal the linear- in-T 
contribution to the low-temperature specific heat and to 
estimate this contribution in order of magnitude, this 
neglection seems permissible (in the limit of small defect 
concentration the cross terms refer to a finite number of 
"optical" vibrations: with the same k but belonging to 
different branches) . This imprecision is the price we have 
to pay for the treatment of a complicated but hermitian 
many-body problem in terms of a set of one-degree-of- 
freedom but non-hermitian subsystems. Anyway, a valu- 
able information has been obtained from this latter ap- 
proach: after solving the mentioned many-body problem, 



the corresponding density of states one must find should 
be similar to that of a damped oscillator [loj ]. 

For further progress, it is necessary to estimate the 
corresponding damping. This estimation, however, is far 
from being trivial even from experimental data. Inelastic 
scattering and nuclear magnetic resonance experiments, 
for instance, permit some estimations. But one must bear 
in mind that these experiments provide a partial infor- 
mation. From the former it can only be extracted the 
high-frequency damping. From the latter it is the low- 
frequency one, but "averaged" over wavevectors. The 
damping of polar optic modes, in particular, can also 
be estimated from experimental data on dielectric losses. 
But these estimates will not be complete: dielectric losses 
are connected to the damping of the long- wavelength op- 
tic modes only. The theory on these dielectric losses is 
well documented, so let us take advantage of this point to 
illustrate that the defect-induced damping does not vanish 
at zero temperature and to make some estimates. 

The problem of the dielectric losses in ferroelectrics due 
to symmetry-breaking defects was theoretically studied 
e.g. in Ref. and reviewed in Ref. 0. The mech- 
anism of loss considered was the defect-permitted radia- 
tion of acoustic waves when applying a time-dependent 
homogeneous electric field (a symmetry-breaking defect, 
e.g. an interstitial atom, induces a local linear coupling 
between the electric field, or polarization, and the strain). 
As a result one finds that the low-frequency damping con- 
stant of a polar mode can be inferred as 
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Here N is the generalized concentration of defects [Tj|, 
u>d is the Debye frequency and n — 0,1,2 for planar, 
linear and point defects respectively. (The defects are 
considered strong defects in the same sense that in Ref. 
[3.) The same results are obtained for non symmetry- 
breaking defects 0, |l| . 

It is worth mentioning that the basic ingredient to 
obtain this defect-induced losses is the accounting for 
the local defect-induced changes in the properties of the 
corresponding crystal, i.e., the local symmetry breaking 
and/or the local changes in the material constants (in- 
homogeneities) . This is sufficient to "connect" the given 
vibrational mode with the acoustic reservoir, which fur- 
ther leads to the corresponding damping (already within 
the approximation of small defect concentration). Hence 
this defect-induced damping is not restricted to polar 
modes: in principle any vibration is affected by these 
defect-induced local changes; and, consequently, Eq. © 
is expected to be valid for non polar modes also. Let 
us stress that the role played by defects, via these local 
changes, is simply to permit the above mentioned con- 
nection between the oscillators and the reservoir (no ad- 
ditional degrees of freedom associated with any "internal 
defect dynamics" are considered). 

As we see in Eq. Q, a frequency-independent damp- 
ing of the "optical" vibrations is obtained due to planar 
defects (n — 0). According to the exposed above, this fur- 
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ther gives a linear-in-T contribution to the specific heat 
which is predominant at low enough temperatures. The 
temperature crossover between this linear-in-T regime 
and the Debye one is T* ~ (Ad 3 ) 1 /^ = (jV p ianard) 1/2 e. 
For a concentration of planar defects corresponding to 
the typical dimensions of the crystal blocks, i.e. one in- 
terblock boundary per fiva. approximately, this crossover 
temperature T* results to be of a few Kelvins. 

The presence of charged defects in ionic crystals pro- 
vide another mechanism of damping for the correspond- 
ing polar modes. The frequency dependence of the di- 
electric losses expected in this case was analyzed in Ref. 
[lCjL Il7j . This analysis can be easily generalized in order 
to compute the frequency dependence of the damping of 
a polar mode with k ^ in the case of a realistic charge 
density distribution. (The choice of the charge density 
distribution in Ref. 0, 0] was somewhat arbitrary, see 
Appendix for further details). As a result, the damping 
coefficient of practically all the polar modes turns out to 
be frequency-independent in the limit u> — > 0. This leads 
to a linear-in-T specific heat at low temperatures, the 
orders of magnitude similar to the obtained above. 

In the case of a frequency-dependent damping such 
that 7 ex U, e.g. for linear defects [n = 1 in Eq. J3J, see 
Refs. [Tj, UK], the damping merely reduces to a renor- 
malization of the optic masses. So in this case, the con- 
tribution of "optical" vibrations will be the "ordinary" 
one. 

Let us finally mention that 7 oc lo 2 can occur due to 
both, point defects [n — 2 in Eq. see Refs. |Llll5|] 
and strongly-correlated charged defects in the case of po- 
lar modes (see Refs. 0, EJ)- Following e.g. Ref. Q 
it can be seen that this further gives a Debye-like con- 
tribution (oc T 3 ) to the low-temperature specific heat 
dgij. Let us stress that this Debye-like contribution has a 
completely different origin than the usual (acoustic) De- 
bye one: it comes from damped "optical" vibrations. It 
is small because of the smallness of the defect concentra- 
tion we are considering. But one can speculate that, for 
high enough concentration of defects, this contribution 
may be comparable to that of the acoustic phonons. 

As an experimental example in which the low- 
temperature contribution to the specific heat due to 
damped "optical" vibrations has probably been observed, 
it is worth mentioning the low-temperature anomaly of 
LiaN reported in Ref. This anomaly cannot be inter- 
preted as a proper glass one (on the basis of the tunneling 
two-level model, for instance) because such a interpreta- 
tion would be inconsistent with the thermal conductivity 
data: the magnitude of this conductivity is too small and 
its temperature dependence is somewhat different from 
that expected in a proper glass (no plateau is observed 
at T > 10 K) j|. But, unfortunately, we are not aware 
of the complementary experiments mentioned above nec- 
essary to confirm (or to reject) that the low-temperature 
specific heat observed in this particular case is due to 
the damping of the corresponding "optical" vibrations. 
Mention that these experiments (inelastic scattering, nu- 



clear magnetic resonance, etc.) must be realized for the 
same sample because the anomaly is expected to be re- 
lated to the presence of defects. Anyway let us notice 
that this low-temperature specific heat anomaly, such 
that C ~ 2 x 1(T 5 J cm- 3 K~ J at T - IK |, can be 
reproduced, in order of magnitude, from Eqs. (J2J) and 
@ with n — 0, a concentration of defects Nd 3 ~ 1CP 4 
(which could correspond to one interblock boundary per 
/im, i.e. N = Api anar d~ 2 ~ 10~ 4 <i~ 3 ), and the typical 
values huJu/kB,fiA/kB ~ 100 K. 

In conclusion, we have shown that the low-temperature 
specific heat of the real crystals may have a significant 
contribution due to damped optical and short-wavelength 
acoustical vibrations. Such a contribution may split into 
linear-in-T and Debye-like terms. The former will prevail 
over the Debye contribution of acoustical vibrations for 
temperatures < 1 K for typical concentrations of defects 
of nominally pure crystals. The latter might be compa- 
rable with the Debye one for high enough concentration 
of defects. Let us stress that these contribution are ex- 
clusively due to damped excitations. This damping is 
due to defects, but defects themselves do not introduce 
any additional degrees of freedom in our considerations. 
If defect excitations are taken into account, some addi- 
tional contribution to the specific heat will be obtained 
similar to that reported for glasses. Mention also that, as 
we have shown, the low-temperature properties of a real 
system may arise from nearly all of its phononic modes 
(all "optical" vibrations of a real crystal contribute to 
the corresponding low-temperature specific heat as we 
have seen) . This does not contradict the general point of 
view according to which, even in glasses, it is accepted 
that only a small number of phononic modes, i.e. the 
long- wavelength acoustic ones, are relevant to the cor- 
responding low-temperature properties. But it indicates 
that the low-energy excitations responsible for non-Debye 
contributions already exist in nominally perfect crystals, 
having their origin in the phononic normal modes. 

We wish to thank S. Vieira, F. Guinea and especially 
M.A. Ramos for useful discussions. 



APPENDIX A: Damping of a polar mode due to a 
disordered charge distribution 

Let us calculate the u> dependence of the damping con- 
stant 7 of a polar mode due to the presence of defects that 
induce a disordered charge distribution in the medium. 
This can be obtained, in the low-frequency regime we 
are interested in, from the u dependence of the imagi- 
nary part of the susceptibility x(w,k): 7(0; — > 0, k) ~ 
UJ ^ 1 x"( UJ ~* 0, k). In the following we shall consider 
that the charge density distribution is correlated in ac- 
cordance with the Debye-Hiickel theory (see, e.g., Ref. 

The results obtained in such a way are different 
from that given in Refs. dill, 

where the choice of the 
charge density distribution was somewhat arbitrary. 
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The k-th Fourier component of the polarization can be written as 

Pi(w,k) = X o(w,k)Bi(w,k) + I ' {plv^iv^e^e^dvdt, (Al) 

where \o is the (scalar) susceptibility in absence of defects (the medium is assumed to be isotropic), V is the volume 
of the system, p is the disordered charge distribution due to the presence of the defects, u is the acoustic displacement 
vector, E is the electric field, and (...) denotes statistical average. The i-th component of this displacement vector 
can be written as [2(J 

tii(w,k) = - Ai - 1 A-. 1 (c,k)^p(k-k')^(c,k'), (A2) 

k' 

where p is the density of the medium and 

A7-V, k) = ^ _ ^2)^2 _ c 2 fc 2) t(^ 2 " C ' fc2 )% + ( C * ~ C t)^] > ( A3 ) 

with q and C( being the longitudinal and transversal velocities of sound respectively. Substituting this expression in 
Eq. l|AT]) we get 

P,(w,k) =xo(a;,k)^( w ,k)- M - 1 ^iir(k-k')Ar. 1 (a;,k')^(a; ) k), (A4) 

k' 

where A'(k) is defined as the spectral density of the charge distribution: 

(p(k)p(k')) oc K (k)5 k ,_ k ,. (A5) 
In accordance with Eq. I|A4|I , the susceptibility due to the defects can be written as 

Xttdrf^.k) = -/i- 1 ^ Jf(k - kOAr.^a;,^). (A6) 

k' 

After integration over wavevectors the non-diagonal terms in X<j',def vanish (we are considering an isotropic medium) . 
So we actually have Xii,def = Xdef<%> where 



Xd. 



rf (w, k) = -m- 1 ^ # (k - k')A^(o;, k') (A7) 



(here double subscript does not imply summation). 

In accordance with the Debye-Hiickel theory (see, e.g., Ref. the charge density around the point r = at 

which the charge density is known, say p(0), is such that 

(p(r + 0)} p(0 ) ~ r- 1 exp(-r/r ), (A8) 

where rjj is the Debye screening length: rjj = [ksT / {Site 2 N)] 1 / 2 , where N is the defect concentration and T can 
be taken as the annealing temperature in our case (the charge distribution is assumed to be in thermal equilibrium 
at this temperature). The desired correlation function (p(O)p(r)) can be obtained by averaging the above magnitude 
over all the possible p(0). This gives 

(p(O)p(r)) = Ne 2 \s(r) - -L- exp(-r/r D )} . (A9) 
L Aitr D r J 

The spectral density of the defect charge distribution is therefore K(k) = K,r 2 D k 2 /(l + rfj/c 2 ), where K, — Ne 2 /V. 

The main contribution to the integral in Eq. I)A7|) then comes from the poles of A^ 1 (w,k'). We are interested in 
the case of small frequencies. Taking into account that the poles of ATT^w, k') are such that k' ~ uj/c (c ~ c/, q), we 
further can distinguish the following cases. For small wavevectors (k <C w/c), in the integrand of Eq. (jA7|l we can 
take 



5 



We then have 

Xdef(u -> 0, k < u/c) oc < 3 , , . _ x (All) 

where wn ~ cdr 1 is the Debye frequency. The u) dependence in these expressions coincides with that reported in 
Refs. [lallS- F° r large wavevectors (k ^> lo/c) , in the integrand of Eq. (| AT|) we can take 

^ (k " k,) ^IT^- (A12) 

We then have 

XdcfV -» 0, fc > w/c) oc Nd 3 (cj/u> D ). (A13) 

From these expressions it follows that, as a result of the presence of charged defects, the damping coefficient of 
practically all the polar modes does not depend on the frequency in the limit to — > 0. 
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